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ABSTRACT 

This paper is one in a series presenting results obtained within the Forma- 
tion and Evolution of Planetary Systems (FEPS) Legacy Science Program on the 
Spitzer Space Telescope. Here we present a study of dust processing and growth 
in seven protoplanetary disks. Our spectra indicate that the circumstellar sili- 
cate dust grains have grown to sizes at least 10 times larger than observed in the 
interstellar medium, and show evidence for a non-negligible (~5% in mass frac- 
tions) contribution from crystalline species. These results are similar to those of 
other studies of protoplanetary disks. In addition, we find a correlation between 
the strength of the amorphous silicate feature and the shape of the spectral en- 
ergy distribution. This latter result is consistent with the growth and subsequent 
gravitational settling of dust grains towards the disk mid-plane. Further, we find 
a change in the relative abundance of the different crystalline species: more en- 
statite relative to forsterite is observed in the inner warm dust population at ~1 
AU, while forsterite dominates in the colder outer regions at ~5 to 15 AU. This 
change in the relative abundances argues for a localized crystallization process 
rather than a radial mixing scenario where crystalline silicates are being trans- 
ported outwards from a single formation region in the hot inner parts of the disk. 
Last, we report the detection of emission from polycyclic aromatic hydrocarbon 
molecules in five out of seven sources. We find a tentative PAH band at 8.2 /xm, 
previously undetected in the spectra of disks around low-mass pre-main-sequence 
stars. 



Subject headings: pre-main-sequence stars 
— planetary systems 



circumstellar disks — circumstellar dust 



1. Introduction 

The circumstellar disks surrounding the 
pre-main-sequence Brown Dwarfs (BD), 
T Tauri stars (TTS) and Herbig Ae/Be 



(HAEBE) systems are believed to be the 
si tes of ongoing pla n et formatio n (e.g . 



Luhman et al. 2007: Natta et al. 2007 



from here on we will refer to these disks 
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as protoplanetary disks). As most of the 
youngest (<1 Myr) s olar-mass sta r s hav e 



circumstellar disks (jStrom et al. 



with typica l masses flBeckwith et al. 



1989h . 



1990) 



and sizes (iMcCaughrean fc O'Delll 11996 
Dutrev et al.lll996h comparable to the 



ex- 
pected values for the primitive solar neb- 
ula, these disks are the natural candi- 
dates for the birth-sites of planets. The 
sub-micron sized dust grains present ini- 
tially in these disks can coagulate to 



like planets (e.g. 


Weidenschilling 




1997; 


Hennine et al. 


2006). Bv deriving the 


com- 



position of the circumstellar dust, and 
identifying the processes governing the 
chemistry and coagulation, valuable in- 
sights can be gained in the workings of 



protoplanetary disks (e.g. iBouwman et al. 



2001, 20081 



van 



Przvgodda et al 



Boekel et al 



2003 



20031. 120051: 



Kessler-Silacci et al 



2006al : ISargent et all 120061 ) and thus the 
planet formation process. The results of 
these analyses can be compared directly 
with the investigation of solar system ob- 
jects like comets, meteorites and interplan- 
etary dust particles (IDPs), which preserve 
a record of the early phases in the evolu- 
tion of the solar system. 

With the IRAS and ISO missions, a 
tremendous advance in our knowledge of 
protoplanetary disks was achieved. How- 
ever, the spectroscopic studies with these 
missions were limited to relatively nearby 
and luminous stars of spectral type A and 
B, and provided only limited knowledge of 
the evolution of solar-mass systems. With 
the launch of the S pitzer Space Telescope 
( Werner et al.l l2004i ). these less luminous 
systems became accessible to observations. 
The FEPS Spitzer Legacy program probes 



the circumstellar dust properties around 
a representative sample of protoplanetary 
disks and debris disks, spanning a wide 
range of circumstellar disk properties and 
covering the major phases o f planet sys 



tern formation and evolution ( Meyer et al. 



20061) 



We present an analysis of the infrared 
(IR) spectra of protoplanetary disks around 
seven pre-main sequence systems observed 
within the FEPS legacy program. These 
are the only systems among the 328 FEPS 
targets showing spectroscopic features 
from solid state dust components. (See 



Meyer et al.l l2006l . for a detailed descrip- 
tion of the FEPS parent sample). In Ta- 
ble [1] the astrophysical parameters of these 
seven stars are compiled. Five of the sys- 
tems, RX J1842.9-3532, RX J1852.3-3700, 
IRXS J132207.2-693812, RX J1111.7- 
7620, and IRXS J161410.6-230542 
identified 



m 



Silverstone et al.l (120061 ) 



were 
as 



optically-thick primordial disks on the ba- 
sis of their IRAC colors and excess emis- 
sion at waveleng t hs sh orter than 8 /xm. 



Silverstone et al.l (120061 ) provide a liter 



ature review for the five stars identified 
in their study. In addition to these five 
systems, HD 143006 and RX J1612.6- 
1859A also show evidence for optically- 
thick disks and solid state features. We 
give a literature review for the latter two 
stars in Appendix A. HD 143006 was 
added to the FEPS sample on the basis 
of its IRAS excess emission, in order to 
search for remnant g as in that system (e.g. 



Pascucci et al.ll2006l ). The solid state fea- 
tures in RX J1612.6-1859A were noticed 
from visual inspection of its spectrum. No 
other sources in the FEPS sample show 
evidence for solid state emission features. 
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Of these 321 sources, most lack excess 
in the IRS wavelengths, and a few ex- 
hibit optically-thin emission from debris. 
No other optically-thick circumstellar disk 
sources are present in the data other than 
the seven with solid state emission features 
discussed here. 

We study the compositional properties 
of the systems listed in Table [1] using 5- 
35 /im spectra obt ained with the Infr ared 
Spectrograph (IRS: lHouck et al.ll2004l ) on- 
board the Spitzer Space Telescope. We fo- 
cus on silicate grain processing (silicates 
make up the bulk of the refractory dust 
mass assuming solar system abundances) 
and perform a quantitative analysis of the 
observed solid state emission. Using a con- 
sistent analysis method for the entire sam- 
ple, the derived grain composition of the 
individual systems can be compared di- 
rectly. This enables us look for signs of 
silicate grain processing. Apart from ana- 
lyzing the 10 /xm spectral region, also ac- 
cessible to ground-based studies, Spitzer 
provides access to the silicate resonances 
at longer wavelengths. This enables us 
to study the dust composition of not only 
the warmest dust (temperatures between 
~ 500 - 1000 K), located in the inner parts 
(~ 0.1 — 1 AU) of the protoplanetary disks 
around solar-like stars, but also the cold 
(~ 100 K) dust component in the outer re- 
gions (~15 AU) of the disks. 

This paper is organized in the follow- 
ing way: In section 2 we discuss the data 
reduction, introduce the main dust com- 
ponents and present the method used to 
analyze the Spitzer IRS spectra. In sec- 
tion 3 modelling results are presented and 
in section 4 we discuss implications for the 
evolution of the dust in the circumstellar 



disks around solar-mass stars. 

2. Observations and Dust Models 

2.1. The IRS Low-Resolution Spec- 
tra 

We obtained low-resolution {R =~ 60 — 
120) spectra with the IRS instrument on- 
board Spitzer. A high accuracy IRS or 
PCRS peak-up (with a la pointing uncer- 
tainty of 0.4" radius) was used to acquire 
targets in the spectrograph slit, thus mini- 
mizing slit losses and assuring high photo- 
metric accuracy to about 10%. Two nod 
positions per cycle were obtained in stan- 
dard staring mode with a minimum of six 
cycles per target for redundancy and to al- 
low the rejection of artifacts introduced by 
bad pixels or cosmic ray hits. The integra- 
tion times were 6s and 14s for the bright- 
est and faintest sources, respectively. The 
targets were observed with the full spec- 
tral coverage of the IRS low-resolution in- 
strument between 5.2 and 35 /zm. Beyond 
35 /xm the spectra suffer from excess noise 
( iHouck et al.ll2004l ) and can not be used. 

Our spectra are based on the droopres 
intermediate data product processed through 
the SSC pipeline S12.0.2. Par tially based 



on th e SMART software package (IHigdon et al. 



2004J . for details on this tool and extraction 
methods), our data is further processed 
using spectral extraction tools developed 
for the FEPS Spitzer science legacy pro- 
gram. As a first step, we correct for the 
background emission and stray-light by 
subtracting the associated pairs of imaged 
spectra of the two nodded positions along 
the slit for each module and order. In this 
way we also correct for pixels having an 
anomalous dark current. Pixels flagged 
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by the SSC data pipeline as "bad" were 
replaced with a value interpolated from 
an 8 pixel perimeter surrounding the er- 
rant pixel. We then extracted the spectra 
from the resulting set of images using a 
6.0 pixel and 5.0 pixel fixed-width aper- 
ture in the spatial dimension for the ob- 
servations with the short- (5.2 - 14 /im) 
and long- wavelength (14-35 /xm) modules, 
respectively. The low-level fringing at 
wavelengths >20 fim was removed using 



theirsf ringe package (ILahuis fc Boogert 
2003h . 



Absolute flux calibration was achieved 
in a self-consistent manner using the en- 
semble of FEPS observations. As the ma- 
jority of the 328 FEPS program stars ex- 
hibit stellar photospheres throughout the 
shorter wavelength range of the IRS, this 
data set is unique for calibrating the IRS 
instrument. From our full set of FEPS ob- 
servations we selected a subset of the spec- 
tra which complied with the following cri- 
teria: 1) colors (24-33 /xm) within la of an 
extra polated best - fit rn o del photosphere 
(e^ iMever et al.l l2004l : ICarpenter et al.l 



20061 ) and 2) Signal-to-noise ratio (SNR) 
larger than 50 with no artifacts within the 
spectra. This resulted for the first order 
of the short wavelength module (SLl) and 
the long wavelength module spectra (LLl, 
LL2) in a set of 16 stars and for the shorter 
wavelength orders (SL2) in a subset of 10 
stars. Together with the stellar photo- 
spheric models for these two subsets, we 
used these spectra to derive the relative 
spectral response functions (RSRFs) for 
the absolute flux calibration. For the rel- 
ative (point to point) calibration, we de- 
rived RSRFs using six calibration stars 
observed at different epochs with the Co- 



hen photospheric models provided by the 
SSC, having a superior SNR. This pro- 
cedure ensures that the RSRFs used to 
calibrate our protoplanetary disk spectra 
have the highest possible SNR and an es- 
timated absolute flux calibration uncer- 
tainty of ~10%. For further details on the 
flux calibration see the explanatory supple- 
ment for the FEPS data, delivered to the 
Spitzer data, archive (Hines et al. 2005). 

The resulting calibrated spectral energy 
distributions are shown in Fig. [H A de- 
tailed view of these spectra in three wave- 
length bands, showing more clearly the 
solid-state emission features is presented 
in Fig. [21 As one can see from Fig. [H a 
wide range in spectral shape and emission 
feature strength can be observed. In the 
next section we will discuss the analysis 
methods and the different silicate compo- 
nents and grain models used to interpret 
the spectra shown in Fig. [2l 

2.2. Outline of Spectral Analysis 
Methods 

For the analysis and interpretation of 
our observations we use a two-fold ap- 
proach. First, we characterize the spec- 
tra by measuring the strengths and po- 
sitions of spectral features. Second, we 
decompose the observed emission features 
to determine the relative contribution of 
different silicate grains. For the spectral 
characterization we employ a method sim- 
ilar to schemes previously used to analyze 
IR spectra fro m the Infrared Spa ce Ob- 
servatory (e.g. iMolster et al.ll2002al]bl') and 

groun d-based observations (e.g. Ivan Boekel et al 
20051 ) . To measure the emission band 
strengths and positions, we flrst flt a low 
order polynomial (second order in the 10 fi 
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region, fifth order for A > 15 /im) to 
the spectral data points excluding obvi- 
ous spectral features (e.g. the continuum). 
As a next step we normalize our spec- 
tra to the fitted continuum as -F;^,norm = 
1 + (i^i.,obs - Fy^cont)/ < Fy^cont >, whcre 
-^i^.obs is the observed Spitzer flux, Fj^cont 
is the fitted continuum, and < -F,y,cont > 
the mean continuum value over the fitted 
wavelength interval. This normalization 
allows us to directly compare the different 
spectra and, ensures that the shape of the 
spectral features remains identical to that 
in the original, unnormalized, spectra. 

To characterize the 10 /im amorphous 
silicate band we determined the fluxes at 
different wavelengths across the silicate 
feature in the normalized spectra. The 
results of this analysis are presented in 
Section 3. Amorphous silicates also show 
a band at around 18 /xm. However, this 
band blends with the (often rising) contin- 
uum emission, as can for instance be seen 
in Fig. [U This makes it difficult, if not im- 
possible, to determine the band strength 
and position with sufficient accuracy for a 
meaningful analysis. We, therefore, opted 
not to use this band in our analysis. At 
longer (>20 /im) wavelengths the emission 
features are mainly produced by crystalline 
silicates. To measure the observed band 
positions and strengths, we simultaneously 
fitted a Gaussian to each of the main spec- 
tral features in the normalized spectra. 
The results of this analysis are also pre- 
sented in Section 3. For the crystalline 
silicate emission bands in the 10 fim wave- 
length region, a slightly modified analy- 
sis method has to be applied. In this 
wavelength range, the emission bands over- 
lap with, and can be dominated by emis- 



sion from small amorphous silicate grains. 
Also, PAH molecules may contribute to the 
observed emission in the short wavelength 
range, overlapping with the emission fea- 
tures of the crystalline silicates. Therefore, 
the crystalline contribution can only be es- 
timated by fitting a detailed grain model 
(discussed in the next section) determining 
the contribution of the amorphous silicates 
and PAH molecules. After subtracting the 
fitted contribution of the amorphous sil- 
icate and PAH molecules, we applied an 
analysis for estimating the different crys- 
talline silicate band strengths in the 10 /xm 
wavelength region, similar to that for the 
longer wavelength crystalline bands. 

2.3. Dust Models 

The grain model and analysis method 
for determining the physical properties of 
the dust grains contributing to the ob- 



served emission are 
cessfuUy applied in 
10 um region (see 


similar to those sue- 
previous studies of the 


Bouwman et al. 


2001; 


van Boekel et al. 


2OO5I. for further discus- 



sion). Here, we also apply this method to 
the silicate emission at longer (~20-30 /im) 
wavelengths. Dust in protoplanetary disks 
has most likely the structure of highly ir- 
regular aggregates containing many differ- 
ent dust constituents, like the interplane- 
tary dust particles collected in the Earth's 
upper atmosphere. Calculating the op- 
tical properties of such complex struc- 
ture has proven to be extremely difficult 



e.g. iHenning Sz Stognienkdll996l ). How 



ever, assuming that the aggregates are ex- 
tremely porous, the individual constituents 
making up the aggregate may interact 
with the radiation field as separate enti 



ties, as is the case of IDPs ( iMolster et al. 
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20031 ). Therefore, we assume that the ob- 
served emission can be represented by the 
sum of the emission of individual dust 
species. Table [2] summarizes the dust 
species used in our analysis, including 



protoplanetary disks (e.er. 


Bouwman et al. 


2001 


van Boekel et al. 


2005; 


Sargent et al. 



20061 ). 

Crystalline silicates such as forsterite 
(Mg2Si04) and enstatite (MgSiOa) have 
many strong and narrow resonances in the 
wavelength range covered by the IRS spec- 
trograph. These rotational/vibrational 
bands of the crystalline dust species al- 
low for an much more accurate determina- 
tion of the chemical composition, grain size 
and shape/porosity compared to the amor- 
phous silicates. We find, like in the above 
mentioned studies, that the pure magne- 
sium end-members of the olivine and py- 
roxene families give the best match to the 
observed spectral features. Based on the 
comparison between the in our spectra ob- 
served band positions and st rengths, and 
laboratory measuremen ts (IFabian et al. 



20011 : iKoike et al.l 120031 ). we find no ev- 
idence for iron containing crystalline sil- 
icates. Still a few open questions re- 
main. Enstatite comes in the form of clino- 
and ortho-enstatite, both having a differ- 
ent crystalline structure. However, both 
forms show very similar emission bands at 
around 10 /im, making it hard to distin- 



guish between the m (e.g. lJager et al.lll998 



Koike et al.l 120001 ). Only at the longer 
wavelengths a clear difference between the 
emission properties can be observed. Un- 
fortunately, as we will discuss in the next 
chapter, no clear enstatite bands can be 
observed at the longer wavelengths, mak- 



ing it difficult to determine the exact en- 
statite structure. We, therefore, have cho- 
sen to use clino-enstatite, allowing for a 
direct compari son between our r e sults and 
the study by Ivan Boekel et al.l (120051 ) of 
a sample of Herbig Ae/Be stars. A sim- 
ilar problem occurs for silica dust grains. 
Silica has nine polymorphs among which 
are quartz, tridimite and cristobalite, dif- 
fering in crystal structure. While of these 
silicates quartz is the most common form 
found on Earth, the common form found 
in interpl anetary dust pa rticles (IDPs) is 
tridimite (lRietmeijerlll988l ) and most likely 
also the most common form in protoplan- 
etary disks. To our knowledge, no sim- 
ilar quality laboratory measurements of 
tridimite exists, enabling us to perform 
an identical analysis as for the other dust 
species. However, if one compares the 
amorphous form of silica and quartz, we 
observes similar bands. It might well be, 
that like in the case of enstatite, it will be 
difficult to tell the different forms of silica 
apart. As quartz is unlikely to exist in pro- 
toplanetary disks we have opted to use an 
amorphous silicate with the stoichiometry 
of silica. 

As alr eady was noted i r i prey ious stud- 



ies (e.g. iBouwman et al.l l200ll ). the ob- 



served emission bands of forsterite, en- 
statite and silica can not be reproduced 
assuming simple homogeneous spherical 
grains, but that one has to adopt an in- 
homogeneous grain structure and/or non- 
spherical grain shape. Here we u se a distri 



butio n of hollow spheres (DHS; iMin et al. 



20051 ) which gives an excellent match to 
the observed band positions and shapes. 
DHS has the advantage over the widely 
used continuous distribution of ellipsoids 
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(CDE; iBohren fc Huffmanlll983r ) that it is 
also defined outside of tfie Rayleigli limit, 
making it possible to investigate the effects 
of grain size on the emission properties of 
the grains. One should realize that the 
above discussed grain models such as DHS 
or CDE should not be taken as an ex- 
act representation of the structure/shape 
of a grain, but rather statistical de- 
scription of the grain shape and structure 
and its deviati on from sphericity and ho- 



mogeneity; (see IVoshchinnikov et al.l 12006 
Min et al. 2006 , for a further discussion of 



grain shape and porosity effects). 

For the amorphous iron/magnesium sil- 
icates the exact composition and grain 
shape can not be constrained as well as 
for the crystalline silicates. Models of 
the broad amorphous silicate resonances 
at ~10 /xm and ~18 /xm are to some ex- 
tent degenerate. Silicate models with vary- 
ing magnesium to iron ratios and differ- 
ent grain shapes/homogeneity can repro 



duce the observed bands (e.g. iMin et al. 



20071 ). To allow for a direct compari- 
son with previous studies, we assume here 
that the amorphous silicates have an equal 
magnesium to iron atomic ratio, and stoi- 
chiometrics consistent with olivine and py- 
roxene. Further, for the amorphous iron- 
magnesium silicates we assume homoge- 
neous, spherical grains for which we calcu- 
late the absorpt ion coefficients using the 
Mie theory (e.g. iBohren fc Huffmanlll983l ) 
and the optical constants listed in Table [2l 

To take into account the effect of grain 
size, we calculated for each of the different 
grain species the opacities for three volume 
equivalent grain radii of 0.1 /im, 1.5 fim 
and 6.0 /im, respectively. These grain sizes 
sample the range of spectroscopically iden- 



tifiable grain sizes at infrared wavelengths 
at the SNR of our data. The resulting cal- 
culated opacities for three successive wave- 
length regimes are shown in Figs. [3^, [3^ 
and [3)d. We present our analysis of the 
various dust species in Sections 13.11 13. 2[ 
and 13.31 

Apart from the silicate dust species, 
PAH molecules also contribute to the ob- 
served emission bands within the spectral 
range of the IRS instrument. These large 
molecules are not in thermal equilibrium 
with the radiation field, but are stochas- 
tically heated. Though many of the ob- 
served spectral features can be linked to 
specific vibrational modes and ionization 
states of the molecules, the exact compo- 
sition remains unclear. As our main focus 
in this paper is the solid-state component 
in the protoplanetary disks, we opted to 
use a template spectrum, based on ob- 
served profiles bylPeeters et al.l (120021 ) and 



van Diedenhoven et al.l (120041 ). for esti- 



mating the contribution of PAH molecules 
to the spectra. We present this analysis in 
Section 13. 4[ 

The thermal infrared emission from pro- 
toplanetary disks is determined by the 
complex spatial distribution of individ- 
ual dust species, and the radial and ver- 
tical optical depth and temperature dis- 
tributions. A completely self-consistent 
model of protoplanetary disks, such as a 
full radiative transfer calculation, would 
also have to account for processes such 
as radial and vertical mixing of circum- 
stellar material and chemical and struc- 
tural alteration of the dust grains un- 
der varying physical conditions in dif- 
ferent parts of the disk. Further, to 
uniquely constrain such models, spatially 
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resolved obs ervations as a functio n of wave- 
length (e.g. iThamm et al.l Il994j ) are nec- 



essary. Instead, we employ a more sim- 
plistic approach which has been demon- 
strated to successfully reproduce the sili- 
cate emission features of HAEBE and TTS 
systems. Detailed self-consistent mod- 



els fe.g. 


Men'shchi 


cov & Hennine 




1997 




Chiane & GoldreichI 


I997I: DuUemond et al. 


2001 


) have shown that the emission fea- 



tures originate from a warm optically thin 
disk surface layer on top an optically thick 
disk interior, producing continuum emis- 
sion. The basic assumption of our simple 
model is that the disk emission features 
can be reproduced by a sum of emissivities 
at a given temperature Tjust, and the con- 
tinuum emission by the emission of a black 
body at a given temperature Tcont, which 
could be interpreted as the typical color 
temperature of the underlying continuum 
at a given wavelength. 

We thus fit the following emission model 
using a linear least square minimization to 
the Spitzer low-resolution spectra: 



— -Bjy (Tcont) Co 



3 5 



I ^ rPAH 



where i?;^ (Tcont) denotes the Planck func- 
tion at the characteristic continuum tem- 
perature Tcont, -Biy(Tdust) the Planck func- 
tion at the characteristic silicate grain tem- 
perature Tdust, i^u^ is the mass absorption 
coefficient for silicate species j (five in to- 
tal) and grain size i (three in total), 1^^^ 
is the PAH template emission spectrum, 
and Co, Cij and Cpah are the weighting 
factors of the continuum, the silicate com- 



ponents and the PAH contribution, respec- 
tively. For the single temperature approx- 
imation to be meaningful, the width of 
the spectral region, and thus the temper- 
ature range, needs to be sufficiently nar- 
row. Therefore, we have split the observed 
spectra into three bands: from 6 to 13 /zm, 
from 17 to 26 /im, and from 26 to 36 /xm. 
Apart from the previous considerations, 
these wavelength boundaries also ensure 
that both continuum and spectral features 
can be well characterized, and that multi- 
ple bands of the main crystalline silicates 
can be modeled simultaneously. This lat- 
ter point is important as it reduces the de- 
generacy problem of overlapping emission 
bands of the different dust species. We 
have limited the maximum temperature in 
our fits to 1500 K, above which all silicate 
dust is expected to be vaporized. 

Fig. [2] shows the resulting best-fit mod- 
els for the three spectral regions. As one 
can see, a good agreement between model 
and observations can be obtained using the 
dust species discussed in this section, with- 
out the need for additional components. 
The resulting model parameters are listed 
in Table [31 We estimate the uncertain- 
ties in the derived band strengths and po- 
sitions and the best-fit dust compositions 
by a monte-ca rlo method as di scussed in 
Section 15.6 of (|Press et al.lll992l ): We gen- 
erate a set of 100 spectra for each obser- 
vation by randomly adding Gaussian noise 
with a 1 (T distribution equal to the error 
in the spectral extraction determined for 
each spectral data point. For each of these 
simulated spectra, we applied the analysis 
methods described above. Our quoted re- 
sults represent the mean values of the de- 
rived quantities and the error is the 1 a 
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standard deviation in the mean. 

3. Results 

Following the dust modelling and spec- 
tral decomposition, our main science re- 
sults are in the areas of dust coagulation 
from the analysis of the 10 /xm amorphous 
silicate feature, dust settling from the ad- 
ditional analysis of the overall SED shape, 
crystallinity from the analysis of the vari- 
ous silicate dust species, and PAHs, which 
remain in the residual emission after sub- 
tracting the identified silicates. 

3.1. Coagulation in protoplanetary 
disks 



Grain size strongly influences the in- 
frared opacities, as can be seen in Fig. [3^ 
to [3t. Typically, for the observations pre- 
sented here, grain sizes (assuming compact 
grains) of up to 10 /im can be inferred 
spectroscopically. This is two orders of 
magnitude larger than the grain sizes typ- 
ically derived for dust grains in the ISM 
(< 0.1 fim). For larger grains, the spectral 
signatures become too weak to distinguish 
them from continuum emission. Amor- 
phous iron-magnesium silicates, making up 
the bulk of the solid state material in the 
ISM, show two broad spectral features at 
about 10 /im and 18 /im from which grain 
sizes can in principle be determined. The 
latter band, however, often blends with 
continuum emission from the disk inte- 
rior and colder dust grains at larger radii, 
making it difficult to determine the grain 
size. For this reason, the 10 /im silicate 
band gives the most accurate information 
concerning amorphous silicate grain sizes. 
Many previous st udies have shown that 
both the shape (iBouwman et al.l l200ll ) 



and strength ( 


van Boekel et al. 


2003 


r 


i005 




Przvffodda et al. 


2003 


: Schee;erer et al 


2006 


Kessler-Silacci et al. '. 


2005. 


2006a) of the 



10 /xm silicate features observed in HAEBE 
and TTS systems are mainly determined 
by the size of the amorphous silicate grains. 

Fig. m shows the correlation between the 
shape and strength of the 10 /xm silicate 
feature as measured by the relation be- 
tween the peak over continuum strength of 
the 10 fim silicate feature and the ratio of 
normalized fluxes at 11.3 to 9.8 fim (filled 
symbols) and at 8.6 to 9.8 /xm (open sym- 
bols). We performed a Kendall's r test on 
the data plotted in Fig. HI This test com- 
putes the probability that (x,y) pairs of 
data are correlated in the sense of the rela- 
tive rank ordering of the x- values compared 
to the relative rank ordering of the y- values 
fjPress et al.lll992f ). A value of Kendall's 



r was computed along with the two-sided 
probability that the variables are not cor- 
related (r = — 1 indicates a perfect anti- 
correlation while r = 1 indicates a perfect 
positive correlation). For the solid points 
in Fig. H Kendall's r = -0.905 with a 
probability P < 0.004 that the data are 
uncorrelated. For the open points in Fig. |U 
we find that r = —0.810 with a probabil- 
ity P < 0.011 that the data are uncorre- 
lated. We interpret these trends as con- 
firmation that the strengths and shapes of 
the silicate features are correlated. Though 
crystalline silicate emission bands will also 
infiuence the shape of the 10 /im silicate 
band, the correlations of Fig. H] can be best 
explained by a change in the typical size 
of the emitting dust grains. The dashed 
line shows the expected correlation for an 
amorphous silicate of olivine stoichiome- 
try with a grain size changing continuously 
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from 0.1 fim (lower right) to 2 fim (upper 
left). The only exception is the 8.6 over 
9.8 /im ratio of RX J1612.6-1859A (#3), 
where the ratio is dominated by an unusu- 
ally strong emission from enstatite and sil- 
ica. 

Fig. [5] shows the relation between the 
mass weighed grain size of the amorphous 
iron-magnesium silicates (open symbols; 
right axis) as derived from our spectral de- 
composition of 0.1, 1.5 and 6 fim grains 
and the peak over continuum value of the 
10 /im silicate band. Examining this rela- 
tion in more detail, we find a Kendall's r 
of -0.81 for the rank ordering of peak-to- 
continuum emission compared to the mean 
size of amorphous silicate grains based on 
the analysis presented in Table [31 The 
probability of these variables not being cor- 
related is 0.011. This correlation demon- 
strates that the peak to continuum ratio is 
correlated to the estimated size of the emit- 
ting amorphous iron-magnesium silicate 
grains. As a note of caution, recall that 
the exact values for the grain sizes derived 
are linked to the grain model used in the 
spectral decomposition. For reasons dis- 
cussed in Section 12.31 we use homogeneous 
spheres to model the amorphous silicates. 
More complex grain structures, like fractal 



ilar results (e.g. 


Voshchinnikov et al. 


2006; 


Min et al. 


2006). but require much larger 



aggregates to reproduce the observed band 
strengths. 

3.2. Coagulation and dust settling 

Having explored grain growth for our 
sample of T Tauri disks compared to the 
ISM grain population, we now ask what 
effect coagulation could have on the disk 



structure itself. As we do not have spa- 
tially resolved observations, we use the 
shape of the SED to infer the geometry 
of the disks. It is long known that the disk 
geometry strongly influences the shape 



of th e SED (e.g. iKenyon fc Hartmann 



19871 ). Strongly flaring disks, intercept- 
ing a substantial fraction of the radia- 
tion from the central star at large radii, 
show a rising SED peaking at around 
100 fim. 'Flattened' disks, on the other 
hand, intercept far less of the radiation 
from the central star, and show a power- 
law-like SED decreasing with wavelength. 



Meeus et al.l (120011 ) showed that the disks 
in HAEBE systems can be divided into 
two gr oups based o n the shape of their 
SEDs. iMeeus et al.1 feoOlf ) interpret this 
bi-modal behaviour of the SED in terms of 
a bi-modal disk geometry, having either a 
flaring or a 'flattened' disk structure. Re- 
cent spatially resolved observations seem 
to confirm the interpretation that the dif- 
ferent SEDs can be linked to d i fferen t disk 
geometries (e.g. iLeinert et al.l 120041 ). As 
an explanation for the different disk ge- 
ometries two different models have been 
proposed: (1) grain growth and gravita- 
tional settling towards the disk mid-plane, 
and (2) 'self shadowing' of the disk surface 
at larger radii by the inner (~ 0.1 AU) 
disk. 

Coagulation models show that if the 
dust grains in the upper layers of a flared 
disk become sufficiently large, they will 
gravitationally settle towards the mid- 
plane of the disk, re sulting in a flattened 
disk geometry fe.g. ISchrapler fc Henning 
2004J : iNomura fc Nakagawall2006h . Recent 



model calculations concerning the effect 
of grain coagulation and settling on the 
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SED show that the changing dust geom- 
etry from a flaring towards a flattened 
geometry due to the grain settling re- 
sults in SEDs consistent with observations 



(D'Alessio et all2006: 


DuUemond & Dominik 


2004b: 


Furlan et al. 


2005, 


2006) 



As an alternative, DuUemond fc Dominik 



(j2004al ) explain the variation in observed 
SED shapes for HAEBE and TTS systems 
as differences in disk heating caused by dif- 
ferences in the self-shadowing of the disk 
surface. If the central star can not directly 
irradiate the outer disk surface due to a 
puffed-up inner disk (causing the surface 
of the outer disk to lay in the shadow of 
the inner disk), the temperature and thus 
the scale height of the shadowed region will 
be substantially lower compared to a disk 
with a directly irradiated surface. This 
lower temperature will result in a smaller 
pressure scale height of the disk and thus 
in a 'flattened' geometry. 

A way to test these models is to look 
for a relation between the measured grain 
sizes and the geometry of the disk. In 
Fig. [6] we explore the relation between 
the amorphous silicate grain sizes, as mea- 
sured by the strength of the 10 fim sil- 
icate feature, and the disk flaring angle, 
as measured by the 30 fim over 13 /xm 
(upper panel) and the 70 fim over 13 /im 
(lower panel) flux ratios. Here, the 13 /xm 
and 30 fim fluxes are synthetic photom- 
etry points derived from our IRS spectra 
using an 1 /xm, respectively, 2 /xm wide box 
centered on the quoted wavelengths. The 
70 /xm fluxes are derived from MIPS obser- 
vations. Fig. M provides quantitative evi- 
dence that these variables are indeed cor- 
related. The Kendall's r for the top panel 
is 0.7 with P = 0.03, while for the bottom 



panel r = 0.81 with P = 0.01, both in- 
dicating a signiflcant positive correlation. 
As we saw in Fig. [5], the peak-to-continuum 
ratio of the 10 /x silicate features is strongly 
anti-correlated with grain size. As a re- 
sult, we conclude from Fig. [6] that grain 
size is strongly correlated with SED slope. 
Larger flaring angles of the disk produce 
larger flux ratios because of enhanced emis- 
sion from cooler grains further out in the 
disk that are exposed to the stellar flux 
by the flaring. With increasing grain size 
(decreasing 10 fim silicate band) the flar- 
ing angle of the disk appears to decrease 
(as traced by the slope of the SED to- 
wards longer wavelengths). Our results are 
in qualitative agreement with the m odel 
predictions of iD'Alessio et al.l (120061 ) and 
DuUemond fc DominikI (j2004bj), and pro- 
vide direct spectroscopic evidence for the 
link between the typical size of the dust 
grains and the disk structure. These re- 
sults suggest that in TTS systems it is co- 
agulation and grain settling towards the 
mid-plane which determines the disk ge- 
ometry and thus the SED, rather then self- 
shadowing by the disk. 

3.3. The crystalline silicates in pro- 
toplanetary disks 

The analysis of features produced by 
crystalline silicates provides important in- 
formation on dust processing in protoplan- 
etary disks. Presently, no crystalline sil- 
icates have been observed in the diffuse 



ISM. iKemper et al.l ( 120051 ) place an up- 
per limit of 1% on the crystalline silicate 
mass fraction of the diffuse ISM. The much 
higher mass fractions observed in proto- 



planetary disks (e.g. 


Bouwman et al.l2001: 


van Boekel et al. 


2005 


, ~5% in this study) 
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therefore, imply that the observed crys- 
talhne sihcates had to be formed in the 
disks themselves, rather then being incor- 
porated directly from the ISM into the 



disk (see lHenning et al.ll2005l ). This makes 



crystalline silicates a tracer of dust process- 
ing and evolution in protoplanetary disks. 
Figs. Hand [H] show the observed crystalline 
silicate bands at mid- and far-IR wave- 
lengths, respectively, for our TTS sample. 
As described in Section 12.21 the spectra 
have been normalized and, for the 10 fim 
regime, corrected for the amorphous sili- 
cate contribution. These normalized spec- 
tra clearly show the emission bands char- 
acteristic of emission from the crystalline 
dust species forsterite, enstatite and silica. 
Using the opacities plotted in Figs. [3K to 
[3b, these three crystalline dust species, in 
combination with the amorphous silicates, 
produce an excellent fit to the Spitzer spec- 
tra as shown in Fig. [2l The overall mass 
fraction of crystalline silicates based on the 
spectral decomposition of the 10 fim sil- 
icate band is around ~5% (see also Ta- 
ble |3]). We find no conclusive evidence 
for other crystalline silicates. Interestingly, 
other crystalline dust species, like iron sul- 
fides or oxides can be expe cted to exist 



in pr otoplanetary disks (e.g. iPollack et al. 



19961 ). Looking at Figs. [2] and [H weak fea- 
tures at 25.5 /im (source number and 5) 
and 30.5 fim (source number and 3) can 
be observed which are not properly repro- 
duced by our dust model. These features, 
however, are 1) not conclusively seen in 
the other targets; 2) the bands do not co- 
incide with known bands of materials like 



iron-sulfide or iron-oxide ( 


Beeemann et al. 


1994; 


Henning et al. 


1995): and 3) no clear 



relation between these weak bands can be 
observed. The possibility that these very 



weak bands are an instrumental artifact in 
the form of some residual fringing can not 
be ruled out at this stage. 

3.3.1. Crystallinity and grain size 

As for amorphous dust grains, the 
shape, strength and wavelength positions 
of the emission bands of the crystalline 
dust species yield information on the typ- 
ical size of the emitting dust grains. As 
the crystalline silicates have many narrow 
spectral features sensitive to shape, size 
and composition of the dust grains, an 
excellent determination of the grain prop- 
erties can be made. Using the analysis 
method described in Section 12.21 Figs. [H] 
and dn] show the correlations between the 
measured position and strengths of the 
main forsterite resonances in our TTs sam- 
ple. No trend in the band positions can be 
observed, which hardly vary from the nom- 
inal position expected for 0.1 /im sized 
grains (filled triangles). This is in sharp 
contrast to the results for the amorphous 
silicate component, where grain sizes of 
at least 6 /im are required to reproduce 
our observations. Also, the band ratios be- 
tween the main forsterite resonances do not 
substantially vary from source to source. 
By comparing the strengths of neighbour- 
ing bands, which minimizes the influence of 
possible source to source variations in the 
mass over temperature distribution, a lin- 
ear relation can be observed. This relation 
reflects differences in the crystalline mass 
fractions in the different sources, rather 
than differences in the grain size of the 
crystalline silicates. A Kendall's r test 
for the data plotted in Fig. [10] results in 
r = 0.81 with P = 0.011 for the 11.3 fim 
versus 10.1 /im band strength, r = 0.90 
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with P = 0.004 for the 33.5 /iin versus 
27.8 fim band strength, and r = 0.24 with 
P = 0.45 for the 23.4 /zm versus 27.8 /im 
band strength. The apparent lack of cor- 
relation between the 23.4 /xm and 27.8 /^m 
band strength most likely reflects the diffi- 
culty in fitting an accurate local continuum 
underneath the 23.4 fim band. Still, those 
sources showing the strongest 23.4 /im 
bands also show the strongest bands at 
27.8 fim. 

From the spectral decomposition in the 
10 /im wavelength region (see also Table[3]), 
we determined the relative mass fraction 
of crystalline silicates contributing to the 
crystalline emission bands. Fig. [5] also 
shows the relation between the observed 
amorphous silicate grain size and the crys- 
talline silicate mass fraction. A Kendall's 
T test gives r = -0.43 with P = 0.18 
suggesting reasonable probability that the 
crystallinity and grain size of the (bulk) 
amorphous dust are uncorrelated. Deter- 
mining the mass fraction of crystalline sil- 
icates depends on how well one can deter- 
mine the mass fraction of the amorphous 
dust component. Within the 10 /xm spec- 
tral region, probing the warmest dust, this 
poses no problem as the amorphous sili- 
cates show a clear spectral signature. At 
the longer (~20 fim) wavelengths, probing 
the colder dust at larger disk radii, how- 
ever, amorphous silicates lack clear spec- 
tral signatures, which makes it difficult to 
determine the physical properties of the 
colder grains. However, the strengths of 
the emission bands over the local contin- 
uum in the ~20 - 30 fim spectral region, 
still provide a measure of the relative crys- 
talline mass fraction function of wave- 
length, and thus as a function of tem- 



perature and hence function of ra- 

dial distance in the disk. Fig. [11] com- 
pares the strength of the 10 /xm silicate 
band -a measure of the typical grain size 
of the amorphous silicates- to the observed 
band strengths and positions of three of the 
main crystalline silicate bands. No con- 
clusive correlation between the amorphous 
and crystalline bands can be observed, im- 
plying that both the observed crystalline 
mass fraction and size of the crystals are 
not correlated with the grain size of the 
amorphous dust, and thus grain growth 
and crystallization appear to be unrelated. 
A Kendall's r test for the relation between 
the peak-over-continuum value versus the 
strength and position of the 27.8 /im band 
results in r = 0.147 with P = 0.65 and 
r = -0.048 with P = 0.88 , respectively. 
For the the relation between the peak- 
over- continuum value versus the strength 
and position of the 33.5 /^m band we find 
T = 0.048 with P = 0.89 and r = 0.52 
with P = 0.1, and for the relation between 
the peak-over-continuum value versus the 
strength and position of the 23.4 fim band 
we find r = 0.33 with P = 0.29 and 
r = 0.62 with P = 0.051, respectively. 
Similar to Fig. [TDl discussed in the pre- 
vious paragraph, the possible weak corre- 
lation we observe for the 23.4 /im band is 
most likely due to a systematic problem in 
correctly determining the local continuum 
using a polynomial fit. As we can see in 
Fig. [H the shape of the SED in the 20 /xm 
region varies considerable. As we shown 
in Section 13. 2[ there is a correlation be- 
tween the shape of the SED and the peak- 
over-continuum value of the 10 /xm silicate 
band. The possible correlation with the 
23.4 /im band, therefore, might simply re- 
fiect a systematic difference in the local 
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continuum estimate between source with 
a towards longer wavelengths rising SED 
and those for which the SED is decreasing. 

At first glance the spectrum of a source 
like IRXS J161410.6-230542 (lower left 
panel in Fig. [2]) seems to exhibit a much 
higher crystalline mass fraction than a 
source like RX J1842.9-3532 (upper left 
panel Fig. [2]). The apparently more pro- 
nounced crystalline silicate bands in the 
former TTS system, however, do not re- 
flect a larger crystalline mass fraction but 
rather less amorphous silicate emission. As 
the amorphous grains grow with respect 
to the crystalline silicates, their relative 
opacities decrease and hence their contri- 
bution to the observed spectra diminishes. 
It is, therefore, crucial to have an accu- 
rate determination of the grain size of the 
(bulk) amorphous silicates. If we would 
have used a smaller maximum grain size in 
our spectral decomposition, we would have 
estimated that no amorphous silicates are 
present in the IRXS J161410.6-230542 sys- 
tem, leading to a 100 % crystalline silicate 
mass fraction. This would, consequently, 
have led to a relation between the grain 
size of the amorphous d ust and the crys- 



tallini ty, as was noted by Ivan Boekel et al. 
(120051 ). This result could be an artifact 



from underestimating the maximum grain 
size in their ground based observations, 
leading to a lower mass estimate of the 
amorphous dust component. Compared 
with these ground-based studies, we should 
note that the Spitzer data have the advan- 
tage to allow for a determination of the 
continuum outside of the spectral window 
accessible from the ground. This results in 
a better characterization of the strength of 
the 10 fim silicate band, and thus of the 



grain size, especially in the case of very 
weak emission features. 

3.3.2. The crystallization process 

Apart from the comparison between 
the crystalline and the amorphous emis- 
sion bands, the inter-comparison between 
the crystalline spectral features can pro- 
vide crucial information. Though we are 
not spatially resolving the protoplanetary 
disks, from model calculations of the tem- 
perature distribution in the disk around a 
typical TTS, one can make an estimate of 
the size of th e emitting region at a given 



wavele ngth (jKessler-Silacci et al.l l2006b 



e.g). iKessler-Silacci et al.l (l2006bl ) argue 



that the silicate feature at around 10 /im 
originates from the disk surface at round 
1 AU from the central star. The spec- 
tral features observed at the longest wave- 
lengths, at around 30 /im, correspond to a 
temperature of ~120 K; the temperature 
at which a black body would emit most 
of his radiation at the given wavelength. 
Assuming a dependency for the dust 
opacity, and a typical stellar temperature 
and radius of 5500 K and 1 Rq^ respec- 
tively, dust grains attain a temperature of 
120 K at a radius of about 15 AU from 
the central star. The implies that with 
the IRS instrument we can probe the in- 
ner ~15 AU of a disk in a TTS system. 
The dust composition of the disks out- 
ward of ~15 AU can not be determined as 
the dust grains at these large radii are to 
cold to emit at wavelengths covered by the 
IRS spectrograph. Comparing the relative 
band strengths of forsterite in the 10 /im 
wavelength region probing the warmest 
dust, with the forsterite bands at longer 
wavelengths probing the colder dust com- 
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ponent further out in the protoplanetary 
disks, a possible trend can be observed (see 
Fig. [T2l) . A Kendall's r test for the rela- 
tion between the 11.3 peak strength 
and the peak strengths of the 23.4, 27.8 
and 33.5 bands results in r = 0.52 
with P = 0.1, r = 0.71 with P = 0.024 
and r = 0.62 with P = 0.051, respectively. 

On average, those systems showing 
the strongest forsterite emission bands in 
the 10 /im spectral region also show the 
strongest bands at the longer wavelengths. 
One could interpret this possible correla- 
tion as evidence of a global crystallization 
process, where the crystallinity of the in- 
ner ~15 AU of the disk is increased from 
non-detectable ISM values to the observed 
mass fractions of about ~5% as measured 
in the 10 /im spectral region. Whether 
this means that the crystalline silicates are 
formed throughout the entire inner disk, 
or formed locally and then redistributed 
throughout the inner disk will be discussed 
in Section HI One should note that due to 
the uncertainties in determining the amor- 
phous mass fractions in a quantatitive way 
at the longer wavelengths (see also Sec- 
tion [2]3]), the radial density profile of the 
crystalline silicates remains uncertain. It 
is possible that in colder parts of the disk, 
at larger radial distances, the fraction of 
crystalline grains deviates from that deter- 
mined for the inner parts. However, our 
results show that, irrespective of the exact 
radial profile, the crystallinity of the very 
inner (~1 AU as measured at 10 /im) and 
outer (~5-15 AU as measured at wave- 
lengths ~ 20 — 30 yum) disk regions are 
related. 

A similar comparison for the enstatite 
emission bands is difficult as at the longer 



wavelengths our spectra show no features 
which can be uniquely identified with en- 
statite. Comparing the strength of the 
main enstatite band at 9.2 /im and the 
longer wavelength bands of forsterite we 
find no clear correlation. This suggests 
that both crystalline species are to a cer- 
tain extent unrelated, and a general in- 
crease of forsterite abundance seems not 
to correspond to a similar increase in en- 
statite abundance. Although the crys- 
talline silicate mass relative to that of the 
amorphous silicates can not be determined 
at the longer wavelengths, our spectral de- 
composition function of wavelength, 
and thus radius of the emitting regions, can 
determine the mass ratios between the dif- 
ferent crystalline silicates. In Fig. [13] the 
relation between the enstatite fraction of 
the crystalline grain populations as mea- 
sured within the 10 /im wavelength region, 
and that determined from the longer wave- 
length spectral regions, is plotted. Sur- 
prisingly, these results imply that the mass 
fraction of enstatite varies in the radial di- 
rection, with a larger enstatite fraction in 
the warmer inner disk than in the colder 
outer (~5 to 15 AU) regions. Another 
difference between the forsterite and en- 
statite grain population is the typical grain 
size. The grain size of the enstatite grains 
(~1 /im) is larger than that of the forsterite 
grains (~0.1 /im). Again this points to- 
wards differences in the formation mech- 
anism and/or conditions producing these 
crystalline species. The implications of the 
above results will be discussed in Section HI 
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3.4. PAH emission in protoplane- 
tary disks around T Tauri stars 

So far, little is known about the pres- 
ence of PAH molecules in disks around 
low-mass pre-main sequence stars. Sys- 
tems with a relatively strong UV field 
can readily excite PAH molecules, as is 
observed for the more massive and lumi- 
nous Herbig Ae/Be systems where clear 
PAH emission signatures can be seen 



(e.g. iMeeus et al.ll200ll : ISloan et al.ll2005f ). 
Though observations of UV poor reflec- 
tion nebula and PAH laboratory mea- 



surements and r no del s (ILi fc Draind 12002 



Mattioda et al.l l2005l ) show that longer 



wavelength photons can also excite PAH 
molecules, little evidence has been found 
for IR emission bands from stochastically 
heated molecules or very small grains from 
the circumstellar disks around low-mass 
stars as studied in this paper. After care- 
fully fitting our silicate model, a residual 
feature centered at 8.2 /im could be seen in 
the model subtracted spectrum of five of 
the observed targets, which we tentatively 
identify with emission from stochastically 
heated PAH molecules. Fig. [H] shows the 
normalized residual spectrum after sub- 
traction of the best-fit silicate model. The 
strongest bands can be seen in the resid- 
ual spectrum of HD 143006 (#2), the most 
luminous source in our sample (G6/8 spec- 
tral type), which also shows clear evidence 
for emission bands at 6.2 fim and 11.2 yum. 

As the PAH molecules require stellar 
photons to be excited, the molecules con- 
tributing to the emission bands have to be 
located at the optically thin surface layers 
of the disk, in direct view of the central 
star. As many papers have s hown for the 
Herbig Ae/Be systems (e.g. iMeeus et al. 



200ll : lAcke fc van den Anckerll2004l ). there 
exists a direct correlation between the 
presence of PAH emission bands and the 
shape of the SED, a measure of the flaring 
angle of the disk. Those disks having an 
SED consistent with a flattened disk struc- 
ture, and hence a less illuminated disk sur- 
face show little or now evidence of PAH 
emission. These results are also confirmed 
by spatially resolved spectroscopy showing 
the PAH emission bands to come from the 
flaring surface at larger radii of the disk 
(Ivan Boekel et al.ll2004bl ). We can observe 
a similar trend, where those sources hav- 
ing SEDs consistent with more flattened 
disk structures show no clear evidence for 
the 8.2 fim feature. This apparent cor- 
relation between the lack of the 8.2 fim 
feature and the flaring angle of the disk, 
therefore, suggests that this emission band 
is not due to an unidentifled dust species in 
thermal equilibrium, but rather a stochas- 
tically heated large molecule or very small 
grain species at the disk surface. 

Surprisingly, the observed band in the 
8 fim region is not centered between 7.7 
to 7.9 fim, where PAH emission is usually 
observed. The band position at 8.2 fim 
has only been me asured in the spectra of a 
few other sources. iPeeters et al.l (120021 ) ob- 
served thi s band position in t wo post-AGB 
stars, and ISloan et al.l (120051 ) observed the 
feature in the HAEBE system HD 135344. 
The observed bands in these three ob- 
jects, however, show a slightly broader pro- 
flle than is observed in the TTS systems. 
The observed relatively weak or absent 
6.2 nm and 11.2 fim features with re- 
spe ct to the 8.2 /xm band is co nsistent with 



the iPeeters et al.l (120021 ) and ISloan et al. 



((20051) studies. Further, the HAEBE sys- 
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tern from the latter study showing the 
8.2 fim PAH band, has the lowest UV-flux 
of the four systems discussed in that pa- 
per, which would be consistent with our 
observations of lower luminosity systems. 
Presently, no conclusive explanation for 
the carrier of the 8.2 /im feature has 
been put forward. We can only speculate 
that the chemistry, ionization state and/or 
structure should be substantially different 
for our sources compared to t he PAH pop- 



ulatio n in HAEBE systems. iLi fc Draine 
(120021 ) argue that in an environment low in 
UV where longer wavelength photons ex- 
cite the PAHs, the PAH molecules have to 
be ionized and/or large to be excited. The 
presumably much stronger X-ray emission 
in the TTS systems with respect to their 
higher mass counterparts may play also a 
role here, driving a different chemistry and 
ionization state. On the one hand. X-ray 
photons may increase the electron abun- 
dance in the gas which would reduce the 
PAH ionization. On the other hand, a 
strong X-ray field could directly (multiple 
times) ionize the PAH molecules due to the 
Auger effect, and thereby change them or 
even destroy the smallest. 



4. Discussion 

The circumstellar disks in T Tauri star 
systems are believed to be the sites of on- 
going planet formation, and thus repre- 
sent an analogue for the proto-solar neb- 
ula. As the disks evolve with time, the 
sub-micron sized dust grains present at 
the formation time of the disks coagu- 
late to form larger objects and eventu- 



ally planet (esimal)s (lBeckwith et al.ll2000 



Henning et al.l |2006| ). By studying the 
characteristics and evolution of the disk 



and its dust composition, valuable insights 
can be obtained into the processes leading 
to the formation of planets, and important 
constraints on disk and planet formation 
models can be derived. Also, by analogy, 
TTS systems can provide clues into the 
early evolution of the solar system. In the 
following we will discuss the implications 
of our findings presented in the previous 
sections. Note that our sample only spans 
a limited range in stellar parameters and 
consists of older, long surviving disks. Our 
conclusions are applicable to similar sys- 
tems, but possible effects of stellar proper- 
ties or age on the evolution of circumstellar 
disks can not be address directly by this 
study. For this our results will have to be 
compared to the disk properties of a larger 
sample, spanning a wide range in stellar 
properties and evolutionary stages. 

4.1. Grain growth and disk struc- 
ture 

Based on our analysis of low-resolution 
infrared spectra of seven TTS systems ob- 
tained with the Spitzer Space Telescope, 
we find clear evidence of processing and 
growth of the sihcate dust species present 
in their protoplanetary disks. We interpret 
observed variations in the thermal emis- 
sion from amorphous silicate dust species 
as evidence for grain sizes which are sub- 
stantially different from those observed 
in the ISM, and argue for grain growth 
within the protoplanetary disks. For the 
first time, we find a clear correlation be- 
tween the strength of the amorphous sil- 
icate feature, measuring the typical grain 
size of these grains, and the shape of the 
mid-infrared spectral energy distribution, 
measuring the disk scale height, i.e., fiar- 
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ing of the disk. In the hterature, two pos- 
sible explanations have been put forward 
to explain the differences in disk flaring: 
First, self- shadowing of the disk due to 
an enhanced scale height of the very inner 
parts. This prevents the central star from 
illuminating the surface of the outer disk, 
required t o sufficientlv heat the disk fo r 
it to flare (jPuUemond &: Dominikll2004al ). 
Second, grain growth and consequent grav- 
itational sett ling toward the disk mid- 
plane (e.g. ISchrapler &: Hennina |2004J: 



Nomur a fc Nakagawall2006l:lD' Alessio et al 



2006 : .DuUemond fc Dominikll2004b : .Furlan et al 
iooil . l2006h . Our results, linking an in- 
crease in grain size to a decrease of the 
disk flaring, clearly argue for the latter 
explanation. 



An earlier study by lApai et al.l (12004 ). 
studying the brown dwarf system CFHT- 
BD-Tau 4, also suggested a correlation be- 
tween the strength of the silicate emis- 
sion band and the sh a pe of the SED. 
Acke fc van den Ancken (120041 ) studied a 
sample of HAEBE systems observed with 
the Infrared Space Observatory (ISO). 
Because no correlation could be found 
between the observed silicate emission 
bands and the shape of the SED, the au- 
thors concluded that their observations 
were consistent with the self-shadowing 
disk model rather than coagulation and 
grain settling. In apparent contradiction, 
Acke et al.l (120041 ) showed a correlation be- 
tween grain size, as measured by the (sub) 
millimeter slope of the SED, and disk flar- 
ing angle, based on the ratio of the near- 
to mid-IR excess. The latter results argue 
for coagulation and grain settling towards 
the mid-plane as the explanation for the 
observed bi-modal distribution of SEDs, 



consistent with our findings here for lower 
mass solar-type stars. 

We suspect that the correlation from IR 
spectroscopy of the derived grain size and 
the disk flaring has been missed so far for 
the HAEBE systems are due to the at least 
one order of magnitude higher stellar lumi- 
nosity of A and B type stars compared to 
the later-type stars studied here. Obser- 
vations of HAEBE systems will, therefore, 
probe a much larger region of the circum- 
stellar disk than the observations of TTS 
or brown dwarf systems for any given wave- 
length, by at least a factor 2 to 3. Given 
that coagulation and settling are a func- 
tion of density and thus radius, gradients 
in the observed grain size and disk struc- 
ture are to be expected. Probing a much 
larger region of the disk, therefore, could 
lead to the loss of any clear correlations In 
the HAEBE systems. 

4.2. The formation of crystalline sil- 
icates 

All systems studied in this paper show 
emission from crystalline silicates, further 
evidence for dust processing within the 
protoplanetary disks as interstellar dust, 
the material present at the formation time 
of the disks, shows no evidence for crys- 
talline silicates. The results of the spectral 
analysis of the crystalline silicate emission 
can be summarized as follows: 

1. All observed targets show emission 
bands of crystalline silicates at both 
higher (~300 - 500 K) temperatures, 
probed in the 10 /im spectral win- 
dow, and lower (~100 K) tempera- 
tures, probed in the 30 fim spectral 
region. 
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2. The observed emission bands are 
consistent with emission from the 
pure magnesium sihcates forsterite 
and enstatite, and sihca. No conclu- 
sive evidence for other species has 
been found. 

3. The average grain size of the crys- 
tals is much smaller than that of the 
amorphous silicates. While amor- 
phous silicates requires grain sizes up 
to 6 /xm, the crystalline emission is 
more consistent with a size of ~0.1 
to 1 fim. 

4. we find no conclusive evidence for a 
correlation between the mass frac- 
tion of the crystalline silicates and 
the grain size of the amorphous sil- 
icates. 

5. The average grain size of the en- 
statite grains as derived from the 
modelling of the 10 spectral re- 
gion is systematically larger than 
that of the forsterite grains. 

6. The strength of the forsterite emis- 
sion features in the 10 fim wave- 
length region and those at the longer 
wavelengths seem to be correlated. 
This could be interpreted as evidence 
that the crystalline mass fraction is 
increased from ISM values to the ob- 
served mass fractions of about ~5%, 
throughout the entire inner ~15 AU 
of the disk. 

7. We find a change in the relative 
abundance of the different crystalline 
silicates from the very inner (~1 AU) 
warm dust population, which is dom- 
inated by enstatite, compared to the 



population at the colder outer disk 
regions at larger radii (~5-15 AU), 
dominated by forsterite. This change 
in relative abundances points to- 
wards a radial dependence of either 
the formation mechanism of the crys- 
talline silicates, or the (non) equilib- 
rium conditions under which they 
formed, and also argues against sub- 
stantial radial mixing of processed 
material from the inner to the outer 
parts of the disk (see paragraphs be- 
low for a detailed discussion). 

How can we interpret the above re- 
sults for the crystalline silicates in a con- 
sistent picture of the physical processes 
occurring in a protoplanetary disk? As 
all our targets show similar crystalline 
mass fractions and compositions, these 
physical processes can be expected to 
occur in all protoplanetary disks. The 
formation of forsterite and enstatite re- 
quires temperatures above ^1000 K (e.g . 
Hallenbeck et al.ll2000l : iFabian et all2000h . 
Such high temperatures are reached near 
the central star where dust grains can 
be readily heated by the stellar radia- 
tion field. In an early active disk phase, 
the accretion energy can provide an ad- 
ditional heat source for the dust in the 
mid-plane regio ns of the inner disk (e.g. 
Bell et al. 2000l ). At these high tempera- 
tures, crystalline silicates can be formed 
either through ann ealing of the amor- 



phous silicates (see IWooden et al.l 12005 



for an overview) or by gas-phase conden- 
sation/annealing a nd gas-solid reactions 



in a cooling gas (IDavis fc Richteii 12003 



Petaev fc WoodI 2005 ). Our observations 
show, however, that a substantial fraction 
of the crystalline silicates have temper- 
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atures (~100 K) far below the required 
crystallization temperature. This implies 
that either large scale radial mixing has to 
occur, transporting the crystalline silicates 
from the hot inner parts to the cooler o uter 
disk (e.g. Bockelee-Morvan et al.l 2002 ). or 
that an additional localised and/or tran- 
sient heating mechanism is operating in 
the cooler regions of the disk, at larger ra- 
dial distances from the central star. Such 
additional heating mechanism might be 
provided by shock-waves in which gas 
and dust grains are heated to the re- 
uired temperatures of above a 1 000 K 
Barker fc Deschll2002l : IPesch et al.ll2005h . 

The question is now if our results can 
distinguish between different formation 
mechanisms and/or locations for the crys- 
talline silicates. The presence of forsterite 
and enstatite, respectively the magnesium 
rich end members of the olivine and pyrox- 
ene solid-solution series, with no evidence 
for crystalline silicates containing substan- 
tial amounts of iron, suggest they formed 
as hig h-temperature, gas-phase conden - 
sates (iDavis fc Richteil l2003l : ICaill l2004h . 
The presence of high-temperature con- 
densates at lower temperatures indicates 
that complete equilibrium condensation is 
not taking place in protoplanetary disks, 
else these species would have been trans- 
formed to other ferromagnesian silicates 
during cooling. Though the exact order 
in which the different dust species con- 
dense out from a cooling gas depends on 
the gas pressure and isolation of the con- 
dens ed dust grains from t he surrounding 
gas (jPetaev &: WoodI l2005l ). most models 
predict forsterite to be condensed out first, 
followed by the formation of enstatite at 
slightly lower temperatures through solid- 



gas reactions between forsterite and Si02 
gas. Iron will condense out slightly before 
or after forsterite as metallic iron grains 
rather than being incorporated into sil- 
icates. The slightly larger grain size of 
enstatite compared to forsterite can be 
the result from the gas-solid reactions and 
the broader temperature range (and thus 
a longer formation window in a slowly 
cooling gas) under which enstatite can be 
formed and is stable. 

An alternative formation mechanism 
for forsterite and enstatite is annealing 
of amorphous silicates without evaporat- 
ing/recondensating the parent grain. The 
annealing process critically depends on 
the structure of the amorphous material 
and the oxygen partial pressure of the 
surrounding gas (i.e. reducing or oxidiz- 
ing conditions). To prevent the forma- 
tion of iron containing olivine and pyrox- 
ene minerals, the amorphous grains have 
to be annealed under reducing conditions. 



Davoisne et al.l (120061 ) show that carbon 
combustion locking up the oxygen in the 
form of CO rather then in FeO, will cause 
the iron the precipitate as metallic parti- 
cles, and not get incorporate d into a sili- 
cate. [Thompson et al.l (120021 ) showed that 
forsterite and silica, independent of the 
stoichiometry of the amorphous material, 
will always form, but that the formation of 
enstatite follows through solid-solid reac- 
tions between forsterite and silica. The 
effectiveness and speed of this reaction 
critically depends on the structure of the 
amorphous grains. In highly porous grains, 
like those formed from smokes, the contact 
surfaces between the individual sub-units 
will be small and, therefore, reactions will 
be slow and require substantial structural 
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modification (like melting) of the amor- 
phous grain. For compact glassy amor- 
phous silicates, the formation of enstatite 
can proceed immediately. Similar as with 
direct condensation, this secondary forma- 
tion of enstatite through forsterite alter- 
ation could also explain why we observe a 
slightly larger grain size for enstatite com- 
pared to forsterite. 

Apart from the above considerations 
two additional observational constraints 
can be imposed to distinguish between 
annealing or condensations as the main 
formation mechanism for crystalline sil- 
icates, namely summary points (3) and 
(4). The implication of these two ob- 
servations is that the simultaneous pres- 
ence of both large amorphous grains and 
small crystalline grains is not caused by a 
preferential transformation of small amor- 
phous grains into crystalline silicates, leav- 
ing the larger amorphous grains. Or with 
other words, the observed change in grain 
size towards larger grains is caused by 
grain growth rather than the removal of 
the smaller amorphous grain population 
trough crystallization in a distribution of 
grain sizes. If the latter would have been 
the correlation between amorphous 

grain size and crystalline mass fraction 
would have been observed. This together 
with the fact that no large crystalline 
grains are observed implies that either the 
crystals form as gas-phase condensates, the 
large amorphous aggregates are disrupted 
before annealing, or that the main crystal- 
lization process occurred before coagula- 
tion. 

Gas-phase condensation can naturally 
explain any difference in grain size between 
the amorphous and crystalline grains. In 



case of an evaporation/condensation zone 
in the disk, where first the amorphous 
grains are evaporated after which crys- 
talhne species can condense out if the 
gas cools, no relation between the orig- 
inal grain size distribution of the amor- 
phous dust and that of the crystalline 
grains will exist. Annealing of large amor- 
phous grains, on the other hand, will lead 
to the formation of large crystals. It is, 
therefore, necessary to disrupt the larger 
aggregates before they get annealed. A 
gradual heating of dust grains as they ap- 
proach the hot inner disk is unlikely to ful- 
fil this requirement. Annealing by shock 
waves could be a more plausible scenario, 
as it could also provide the mechanism to 
disrupt the larger amorphous aggregates. 
This formation mechanism would be simi- 
lar as the formation scenario of chondrules 



e.g. iDesch et al.ll2005f ) 



In case of an early formation time of 
the crystalline silicates, that is before co- 
agulation and during the high accretion 
phases, viscous dissipation can efficiently 
heat the disk u p to several AU from the 
central star (e.g Bell et al. 20001 ). An early 



formation time would imply that the crys- 
talline silicate content of protoplanetary 
disks would remain fixed during their fur- 
ther evolution. This is consistent with 
previous studies who reported no correla- 



and systemic ages (e.g. 


van Boekel et al. 


2005 




Kessler-Silacci et al. 




2006a. 


2005 





though any correlation might also be lost 
due to the large uncertainties in age de- 
terminations. On the other hand, all sys- 
tems studied in this paper show crystalline 
silicates. Studies of on average younger 
samples show that a substantial fraction 



22 



of dis ks show no crystallin e silica tes (e.g. 
50 % iKessler-Silacci et all bnOBah . This 
would argue in favour of crystaUinity in- 
creasing with time. Further, recent studies 
of FU Orionis objects show no sign of crys- 
talline silicates, independent of whether 



envelopes ( 


Ouanz et al. 


2006; 


Green et al. 


2006 




Quanz et al. 20071). If the notion 



is correct that the FU Orionis type ob- 
jects represent the early high accretion 
stages of TTS evolution, crystalline sili- 
cates should be present if the bulk of the 
crystalline silicates would form during the 
early evoluti o n of t he disks. Interestingly, 
Quanz et al.l (120061 ) show that FU Orionis 
already shows signs of grain growth. If co- 
agulation sets in this early, an alternative 
scenario for the formation of crystalline 
silicates is required. 

Note that there is no difference in co- 
agulation behaviour between crystalline 
and amorphous silicates. The difference 
in grain size between the crystalline and 
amorphous silicates can be interpreted in 
terms of the optical properties of larger 
composite aggregates. Suppose a popula- 
tion of small, predominantly amorphous 
dust grains coagulates into larger aggre- 
gates. These larger aggregates will obvi- 
ously consists mainly of amorphous ma- 
terial but with few, from each other iso- 
lated crystalline grains which will there- 
fore react with the radiation field as sep- 
arate entities. The optical properties of 
such composite aggregates will likely re- 
semble that of the combined properties 
of a larger amorphous grain with smaller 
crystalline particles, similar as is observed 
in cometary spectra and IDPs. 



The two remaining summary points, (6) 
and most importantly (7), to be discussed 
here, place further constraints on the for- 
mation mechanism of the crystalline sili- 
cates and on the location in the disk where 
the crystalline silicates are forming. On 
the one hand our data suggests a general 
increase in the crystalline mass fraction in 
the inner ~15 AU of the disks, on the other 
hand we observe a gradient in the crys- 
talline dust composition within this inner 
region. The changing forsterite over en- 
statite mass ratio points to differences in 
the condensation or annealing conditions. 
As discussed in the previous paragraphs, 
the formation of enstatite follows after the 
formation of forsterite. The absence or 
much lower enstatite mass fraction derived 
from the longer wavelength bands, prob- 
ing the cooler dust at about 5-15 AU com- 
pared to the mass ratio derived from the 
10 /im wavelength region, shows that much 
less or even none of the forsterite is con- 
verted into enstatite at larger radial dis- 
tances from the central star. In the case 
of gas-phase condensation as the formation 
mechanism for the crystalline silicates, this 
implies non-equilibrium condensation con- 
ditions at larger disk radii, in contrast to 
the very inner disk regions. If forsterite 
forms locally at ~5-15 AU in the disk 
it has to be isolated from the surround- 
ing SiO gas before they react to form en- 
statite. This could be achieved by coagu- 
lation, or by rapid cooling of the surround- 
ing gas, both inhibiting the gas-sohd reac- 
tions leading to the formation of enstatite. 
Alternatively, in case forsterite is formed 
through annealing of ISM material, the low 
enstatite to forsterite ratio implies a rather 
porous structure of the amorphous silicates 
and brief heating events, like shocks, pre- 
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venting the formation of enstatite. 

The high temperature region in the in- 
ner disk near the central star, is the natu- 
ral location where crystallization must oc- 
cur. Assuming that the inner disk region 
is the sole region where crystalline silicates 
are formed, the presence of crystalline sili- 
cates at low temperatures implies that ef- 
ficient radial mixing of the processed dust 
from the inner disk to larger radii must oc- 
cur. At the early, high accretion phases of 
the disk, the transport of angular momen- 
tum is expected to efficiently mix mate- 
rial from the hot inner to the cold outer 



parts of the disk ( 


Bockelee-Morvan et al. 


2002; 


DuUemond et al. 


I2OO6 


). Though at 



a first glance this scenario would be able 
to explain our observations, a strong argu- 
ment against it is the radial dependency 
of the crystalline silicate composition: If 
the crystalline silicates originates from a 
single formation region before being dis- 
tributed through the disk, the composition 
of the processed material would be con- 
stant with radius. Indeed, detailed disk 



models by lGaill (120041 ). taking into account 
radial mixing, predicting the dust compo- 
sition as a function of radius, s how e xactly 
this counter argument. While iGaill (120041 ) 
predict the composition of the crystalline 
silicates in inner regions of the disk to 
be dominated by enstatite, consistent with 
our observations, the radial mixing model 
will also predict a large mass fraction of en- 
statite at larger radii, in contradiction with 
our results. The only way to save the radial 
mixing scenario is if during the active mix- 
ing phase, non-equilibrium conditions pre- 
vailed in the inner parts of the disk, form- 
ing only forsterite. This might be achieved 
by very rapid outward transportation of 



the formed forsterite preventing prolonged 
heating of the grains or substantial gas- 
solid reactions which could lead to the for- 
mation of enstatite. During the passive, 
low accretion disk phases, the phase were 
we are observing our TTS sample, were ra- 
dial mixing is not expected to be efficient, 
equilibrium conditions might be reached 
again, resulting in the observed large frac- 
tion of enstatite in the inner disks. At 
present, however, there is no theoretical 
foundation for this scenario. 

The alternative model to the radial mix- 
ing scenario is the in-situ formation of 
crystalline silicates at radii of about 5- 
15 AU. As at these distances stellar ra- 
diation can not heat the dust to suffi- 
ciently high temperatures for crystalliza- 
tion to occur, an alternative heating mech- 
anism is required. A promising model is 
the loc al heating of dust and gas by shock 



waves ( Harker fc DeschI l2002l : iDesch et al. 



2OO5I ). For typical densities and shock 
speeds, small dust grains should be ef- 
ficiently heated to high temperatures up 
to disk radii of about 10 AU. A possible 
source for the shocks could be gravitational 
instabilities or planetary mass companions 
within the disk. Shock heating could pro- 
vide sufficiently brief heating events such 
that equilibrium conditions or prolonged 
annealing will not occur, preventing the 
formation of enstatite at larger radii. In 
the very inner (~1 AU) disk, at higher 
densities, apart from eating by the central 
star, cooling can be expected to be longer 
such that equilibrium conditions could be 
reached. The suggested overall increase 
of the forsterite mass fraction in the in- 
ner ~15 AU of the disks, might reflect 
that spiral waves could affect the entire 
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disk region, causing shock to occur over 
a large range of radii or simply that the 
forsterite, produced at 5-15 AU from the 
central star, is accreted inwards to smaller 
radii. This shock-wave scenario could ex- 
plain why some systems do not show any 
evidence for the presence of crystalline sil- 
icates at th e longest or even a l l wave - 
lengths (e.g. Kessler-Silacci et"aI1 2006a ): 
these systems might not (yet) have pro- 
duced disk instabilities or sufficiently large 
planetary sized objects to produce shock- 
waves. 

Ideally, the spatial distribution of the 
different silicate species, should be di- 
rectly determined by spatially resolved 
spectroscopy. Due to observational dif- 
ficulties, this is at present limited to 
the much brighter HAEBE and the most 
nearby TTS systems. Using mid-IR inter- 
ferometric observations obtained with the 



VLTI/M IDI instrument, Ivan Boekel et al. 



(l2004al ) could spatially separate contri- 
bution to the IR emission coming from 
the inner ~2 AU of the disk, and that 
coming from larger radii, for a sample 
of HAEBE systems. These observations 
clearly showed that the crystalline sili- 
cates in the observed HAEBE systems 
are mainly concentrated in the inner, high 
temperature, parts of the disk, and have 
a composition consistent with a forma- 
tion under (near) equilibrium conditions. 
However, before concluding that all crys- 
talline silicates form in the inner parts of 
the disk, one has to real ize that of the by 
van Boekel et al.l (l2004al ) studied systems, 
show no or substantially weaker crystalline 
silicate bands at the longer (20-30 fim) 
wav elengths in their Sp itzer and ISO spec- 
tra (IMeeus et al.ll200ll . Bouwman et al, in 



prep.). This suggest that these HAEBE 
systems, in contrast to the in disks in our 
TTS sample, did not witnessed substantial 
radial mixing or in-situ formation of crys- 
talline silicates at larger radii, with only 
the very inner regions of the HAEBE disks 
containing observable quantities of crys- 
talline silicates. 

4.3. PAH molecules in disks sur- 
rounding T Tauri stars. 

Though the main focus of this pa- 
per is the silicate dust processing we re- 
port the tentative detection of emission 
bands from polycyclic aromatic hydrocar- 
bon molecules. Emission bands of PAH 
molecules are commonly found in envi- 
ronments with a strong UV flux like the 
circumstellar disks around HAEBE stars. 
As PAH molecules can contain a substan- 
tial fraction of the available carbon, and 
are spectroscopically easier detectable as 
carbon in larger grains, they can provide 
us with important clues concerning the 
carbon chemistry of circumstellar environ- 
ments. Characterizing the PAH popula- 
tion is also important for understanding 
the observed gas temperatures in disks, 
as the stochastically heated molecules 
strongly influence the temperature through 
photo-electric heating. Since the gas tem- 
perature determines the pressure scale- 
height of the disk, the presence of PAH 
molecules can, therefore, have a substan- 
tial influence on the disk geometry. Five 
out of the seven TTS systems of our sam- 
ple show a band at 8.2 /xm we identify with 
emission from PAH molecules. This is the 
first time this band is observed in low- 
mass pre-main-sequence systems. The rel- 
ative high fraction of systems in our sample 
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showing evidence for PAH molecules seems 
to be contradicting other studies such as 
the Core to Di s ks Sp itzer legacy study 



by (iGeers et al.l l2006l ). which found no 



PAH emission around sources with spectral 
types later than G8. However, many stud- 
ies are based on the analysis of the 6.2 and 
11. 3 jUm PAH ba r ids. A s already noticed 
by IPeeters et al.l ((2002), these latter two 
bands seem to be suppressed in the PAH 
population producing the 8.2 /im feature. 
Just based on the 6.2 and 11.3 /xm features, 
we would only have detected PAH emission 
in the HD 143006 system, a G6/8-type 
star, which wou ld be consistent w ith the 
conclusions by (iGeers et al.l l2006l ). Also, 
as the 8.2 fim feature coincides with the 
strong amorphous silicate band, one can 
only make a firm conclusion concerning 
its presence after carefully modelling the 
silicate emission bands. The fact that we 
observe a PAH band at 8.2 /im rather than 
the commonly observed band position be- 
tween 7.7 to 7.9yum suggests a fundamental 
difference in carbon chemistry between our 
sample and the more luminous HAEBE 
and F- and early G-type TTS systems. 
At this point we can only speculate that 
the difference in carbon chemistry could 
be linked to differences in X-ray fluxes or 
stellar wind properties of the intermediate- 
mass systems compared to those of low- 
mass pre-main-sequence stars. If so this 
could also influence the possibility and 
chemical path by which more complex pre- 
biotic organic molecules can be formed in 
protoplanetary disks. 
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A. HD 143006 



HP 143006 is a w e ll-studied fe.g. 



Sylvester et all Il996l : ICoulson et al. 



Garcia-Lario et al. 


1990l: 


Walker & Wolstencroft 


1988 


1998: Natta et al. 


2004J: 


Dent et al. 


20051) classical 1 



Tauri star of spectral type G6/8 which appears to be member of the Upper Scorpius OB as- 
sociation (jMamajek et al.ll2004). It d isplays many T Tauri c haracteristics inclu ding infrared 
to millimeter excess (IOdenwaldlll986[). optical emission lines ( ISt ephensoni 1 1 9 86l ) . and coronal 



activity revealed through x-rays (ISciortino et al 
system, with an age of 3-5 Myr. 



19981 ). It is thus an "old" accretion disk 



B. RX J1612.6-1859A 



RX J1612.6- 1 859A and RX J1612.6-1859B are the names assigned in the nomenclature 
of iMartin et all (119981 ) for a pair of emission line sources in the p Oph vicinity. The former 
object is an MO star also known as 2MASS J16123916-1859284 or GSC 6209-01312 that has 
been studied only by optical spectroscopy (IMartin et al.lll998l : lThelll964j). The latter obj ect, 
by contrast, is a well-kno wn T Tauri star, also known as ScoPMS 52 (j Walter et all 1 19941 ) or 
Wa Oph 3 (jWalteiill986l ). This KO star was, in fact, the target of the FEPS observations 
and the source discussed in our companion paper (ISilverstone et al.ll2006l ) as not having any 
excess emission out to 8 yum (and indeed not out to 24 /xm; Carpenter et al., in preparation). 
For our IRS observations of this source, however, the IRS peak-up routine centered on to 
the mid-IR bright object RX J1612.6-1859A in the slit rather than RX J1612.6-1859B 20" 
away. The observations discussed herein are thus of RX J1612.6-1859A, an MO emission line 
object. 
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Table 1 

astrophysical parameters of program stars. 



ID# 


Name 


Oi 


5 


d 


log(Age) 


Sp. Type 


Ref. 




Av 








(2000) 


(2000) 


[pc] 


[y] 




[K] 




[mag] 


[L0] 





RX J1842.9-3532 


18:42:57.98 


-35:32:42.73 


145 


6.63 (4 Myr) 


K2 


2,3,1 


4995 


1.03 


1.0 


1 


RX J1852.3-3700 


18:52:17.30 


-37:00:11.93 


145 


6.5-7 (4 Myr) 


K3 


2,3,1 


4759 


0.92 


0.6 


2 


HD 143006 


15:58:36.92 


-22:57:15.35 


145 


6.70 (5 Myr) 


G6/8 


7,8,9 


5884 


1.63 


2.5 


3 


RX J1612.6-1859A 

(2MASS J16123916-1859284) 


16:12:39.18 


-18:59:28.0 


145 


6.70 (3 Myr) 


MO 


11 


3800 


1.80 


0.5 


4 


IRXS J132207.2-693812 
(PDS 66) 


13:22:07.53 


-69:38:12.18 


86 


7.23 (17 Myr) 


KllVe 


10 


5228 


1.22 


1.3 


5 


RX Jllll.7-7620 


11:11:46.32 


-76:20:09.21 


163 


6.69 (5 Myr) 


Kl 


5,3,4 


4621 


1.30 


1.6 


6 


IRXS J161410.6-230542 


16:14:11.08 


-23:05:36.26 


145 


6.70 (5 Myr) 


KO 


7,8,6 


4963 


1.48 


3.2 



([PZ99] J161411.0-230536) 



References for distance, age, spectral type — (1) Neuhauser et al. 2000; (2) Neuhauser & Forbrich 2007; (3) Hillenbrand et al. 2007, in 
preparation; (4) Alcala et al. 1995; (5) Luhman 2007; (6) Preibisch et al. 1998; (7) de Zeeuw et al. 1999; (8) Preibisch et al. 2002; (9) 
Houck & Smith-Moore 1988; (10) Mamajck ct al. 2002: (11) Martin et al. 1998. 

Notes for temperature, extinction, luminosity - Spectral types are from optical spectroscopy as cited above. B-V and V-K colors are used 
in conjunction with spectral types to estimate effective temperatures and visual extinction values; see Carpenter et al. 2007, in preparation, 
for details. In the case of RX J1612.6-1859A J-H and H-K colors are used and luminosity is computed from J-band bolometric correction. 



Table 2 

Overview of dust species used. For each component we specify its lattice 
structure, chemical composition, shape and reference to the laboratory 
measurements of the optical constants. for the homogeneous spheres we 
used mie theory to calculate the opacities. for the inhomogeneous 



SPHE RES, W E USED THE DISTRIBUTION OF HOLLOW SPHERES (DHS: IMIN ET AL 
20051 ). TO SIMULATE GRAIN DEVIATING FROM PERFECT SYMMETRY. 



# 


Species 


state 


Chemical 
formula 


Shape 


Ref 


1 


Amorphous silicate 


A 


MgFeSi04 


Homogeneous 


(1) 




(Olivine stoichiometry) 










2 


Amorphous silicate 


A 


MgFeSizOe 


Homogeneous 


(1) 




(Pyroxene stoichiometry) 










3 


Forsterite 


C 


Mg2Si04 


Inhomogeneous 


(2) 


4 


Clino Enstatite 


C 


MgSiOs 


Inhomogeneous 


(3) 


5 


Silica 


A 


Si02 


Inhomogeneous 


(4) 



References. — ( i jlDorschner et all hm^y f2ls ervoin &: Pirioul 
rs ljager et all (|l998l ): V4 iHenning fc Mutschkel (|l997l i: 
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Table 3 

The best fit values of the parameters in our compositional fits to the 8 
to 13 fxm spectral region. the abundances of the various dust species are 

GIVEN AS A PERCENTAGE OF THE TOTAL DUST MASS, excluding THE DUST RESPONSIBLE 
FOR THE CONTINUUM EMISSION. AlSO LISTED IS THE MASS AVERAGED GRAIN SIZE OF 
THE DIFFERENT DUST SPECIES. IF A SPECIES WAS NOT FOUND, OR UNCONSTRAINED BY 

the spectra, this is indicated by a - symbol. the pah and continuum flux 
contributions are listed as percentages of the total integrated flux over 
the 10 jlm region, contained in these components. these are measures for 
the relative flux contributions, but cannot be interpreted as relative 
dust masses. the resulting best fit model spectra are plotted in fig. [5] 

(light grey lines). 



# 






Td^st [K] 


Cont. 


Amor ph. Olivine 


Amor ph. Pyroxene 


Forstcrite 


Enstatitc 




Silica 


PAH 










contr. 


mass < a > 


mass < a > 


mass < a > 


mass < a > 


mas 


is < a > 


contr. 



8-13^m spectral region 



11.1±0.5 1500. Oto 5! 

1 3.4±0.3 1500. otoS 



.+3.9 



2 17.2±0.7 1459. 0"' 



257.4^^ . 
296.2^2:6 



n+17.1 

-13.5 

3 8.8±0.7 1030.6^94 7 369.9155'* 
a+2.7 ■ 
-14.6 
n+0.0 
-0.0 



4 11.1±0.5 778.8 

5 4.0±0.4 1500. 0;! 

6 6.9±0.4 617.0^7^1^ 



44.0 
183.81^4=8 
236.6tl^| 
305.1"' 



1.7-2B[im spectral region 

iis.o 

1 7.0±0.7 105.3 



1+5.6 
-5.8 



0.5214+[5- 
0.6361^0 
0.6255+[5' 
0.6792+°' 
0.6054+°' 
0.8544+°' 
0.8206+°' 



0009 
0009 
0027 
0025 
0011 
0013 
0263 
0210 
0031 
0011 
0011 
0010 
0055 
0096 



0.7902^0; 
0.9872^°; 
0.4225t°; 
0.2933^°; 
0. 63691°; 
0.8236t°; 
0.5862^°; 



0063 
0063 
0042 
0058 
0120 
0053 
2511 
2375 
0367 
0412 
0138 
0160 
0810 
1153 



1.16' 
3. 

0.13 
6.00 
5.57' 
6.0 
6.0 



+ 0.09 
0.09 
+ 0.07 
-0.13 
+ 0.47 
0.03 
+ 0.0 

0.009 
+ 0.08 
0.05 
+ 0.0 
0.0 
+ 0.0 
0.0 



0.1787^0 
0.0016^°' 
0.5296^°' 
0.6452+°' 
0.3279^°' 
0.1060+°' 
0.3728+°' 



0071 
0072 
0064 
0014 
0047 
0137 
2366 
2529 
0404 
0370 
0153 
0133 
1200 
0802 



,,+0.09 
-0.08 
1+0.0 
-0.0 

n+0.06 

-0.05 
5+0.34 
-1.94 
7+0.45 
-0.75 
1+0.0 
-0.0 



1.12 
0.1 
5.30 
4.85 
2.27 
0.1 

5 999+° °°l 



0. 0146^2; 

0. 00161°; 

0.0149t°; 

0.0141^2; 
0.0114^°; 

0.0272^°; 
0.0058^2; 



0004 
0004 
0011 
0009 
0003 
0003 
0044 
0018 
0003 
0005 
0025 
0023 
0005 
0005 



+0.0 
0.0 
+0.0 
0.0 

1+°° 
'J'^-0.0 
+0 



o.i-" 
ai-i 



0.15 
0.1 
0.97 
0.1 



+0.0 
-0.0 
+0 


+0.0 
-0.0 



0.0165^0 
0. 0091^0 
0.0330l° 
^ni 0.0342t° 
0.0238^0 
ol 0.03491° 
0.0249+° 



0016 
0015 
0036 
0038 
0010 
0011 
0061 
0054 
0013 
0026 
0023 
0022 
0022 
0024 



1.41 
1.496 
1.13 
0.97 
1.39 
1.21 



+0.08 
-0.09 
+0.08 
0.32 
+0.004 
0.02 
+0.13 
0.13 
+0.04 
0.07 
+0.07 
0.08 
+0.08 
0.09 



0006 + °'°°^° 



0133+°'°°*° 
u.uiaa_o (,(,27 



n nn»-?+° °°i9 
0.0083_o (,(,^5 

0104+°'°°°^ 
"•"-^"^-0.0009 



c+O.O 
^-0.0 



0.0066t°„;ro 
0.0082j:°;™l 
- + 0.000 
-0.000 



0.0246^ 



37+0-27 (J (j,(j3+0.003 
'^"^'-0 22 ^■^-'-^■^-0.002 
0126+°'°°° 
35+°'=^ 0022+°'°°° 

1+°'° 
"■^-0.0 



0066+°'°°° 
"•""""-o.ooi 



28.3±1.4 46.3ti"4 

,+30.1 
-0.6 

2 18.8±1.0 33.0t|'| 
,+6.0 
-3.2 



3 5.9±0.6 18.2"f 



4 13.0±0.8 106.0^74 

5 5.3±0.5 177.0t2;6 



6 3.9±0.5 340.0 



+ 9.7 
110.6 



+0.0 
-0.0 
84.95tg;°= 
112.7j:°'| 
197.7l^-°4 
152.61=3'^ 

294.5t^3t;^3 
+ 71.6 
13.1 



155.1"' 



0.0321+°' 
0.2529^°' 
0.0096^°' 
0.00131°' 
0.1748^0 
0.7317+°' 

,+o: 



0.6692^ 



0022 
0018 
0186 
0289 
0047 
0024 
0010 
0006 
0513 
0388 
0120 
0104 
0308 
1003 



0.9592"f 



+0.0066 
0.0071 
0.8650t°„;°„lf, 
0.7384l°;Jf,° 

l214+° °348 
°'^^^*-0.0353 
n ■?a/iK+°-°°12 
0.3846_(, 0554 



i.s: 

1.61"* 



+ 0.0 
-0.0 
+ 0.14 
0.09 
1 5+°'° 

-'-■^-o.o 

1 5 + °'° 
1 5 + °'° 



0.0328t°„ °^5t ^■Ot° ° 



0.0 



0.016ll°' 
0.1141+°' 
0.2459^°' 
0.8628^°' 
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Fig. 1. — The spectral energy distributions the the TTS systems observed within the FEPS 
legacy program. Shown are the Spitzer low-resolution spectra scaled to the adopted distance 
and luminosity of the individual stars as listed in Table [H For clarity the spectra are off-set 
from each other, ordered from top to bottom based on the observed slope of the SED, by 
multiplying by 20, 5, 0.75, 1, 2.2, 0.25 and 0.2, respectively. 
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Fig. 2. — Spitzer low-resolution spectra of 7 TTS observed within the FEPS legacy program. 
The left panels show the spectra in the wavelength region between 8 and 13 fim. The central 
panels show the spectral region between 17 and 26 fim, and the right panels the spectral 
region between 26 and 36 fim. The spectra are ordered from top to bottom according to the 
decreasing strength above continuum of the 10 /im silicate band. Also plotted in this figure 
is a compositional fit to the spectra (light grey lines). For details on the model fits see 12.31 
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Fig. 3. — To unity scaled opacities for the 
different dust species listed in Table El Fig- 
ure 3a shows the opacities between 8 and 
13 iim, Figure 3b between 17 and 26 /xm, 
and Figure 3c between 26 and 36 fim. In 
each panel the solid lines represent the 
opacities for a grain size of 0.1 fim, the 
dashed lines the opcities for a 1.5 /xm grain 
size, and the dot-dash lines the opcities of 
6 yum sized grains. Also indicated in each 
panel are the scaling factors for each of the 
opacity curves, used in the normalization 
to unity. 
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Fig. 3. — Continued. 
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Fig. 4. — Correlations between the strength and the shape of the 10 //m sihcate band. The 
filled symbols show the correlation between the peak over continuum ratio of the 10 fim 
silicate band and the ratio of the normalized flux at 11.3 fim over 9.8 /im (left axis). The 
open symbols show the correlation between the peak over continuum ratio and the ratio of 
the normalized flux at 8.6 /im over 9.8 ^m (right axis). Note that the formal error on the 
plotted quantities is smaller then the size of the used symbols. The numbers correspond to 
the ID numbers of our target stars as listed in Table 1. Also plotted in this flgure, represented 
by the dashed lines, is the calculated behavior for both flux ratios of amorphous olivine grains 
for a continuous changing grain size from 0.1 /xm (lower right) to 2.0 /i,m (upper left). 
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Fig. 5. — Correlation between the derived mass fraction of crystalline silicates (filled symbols; 
left axis) and the peak-ovcr-continuum ratio of the 10 /im silicate band and the mass averaged 
grain size of the amorphous silicates (open symbols; right axis) emitting at this wavelength 
region. As one can see, as the amorphous grains become bigger the 10 fim silicate band 
becomes weaker. No correlation between the mass fraction of crystalline silicates and the 
typical grain size of the amorphous silicates can be observed. 
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Fig. 6. — Correlation between the far-IR slope of the spectral energy distribution, a measure 
for disk geometry, and the peak-over-continuum ratio of the 10 /im silicate band, a measure 
for the typical grain size. The far-IR slope is measured by the ratio of the flux at 30 /im over 
13 (upper panel) and the 70 fim over 13 fim flux ratio (lower panel). A clear correlation 
can be seen: As the 10 /im silicate band becomes weaker, the slope of the spectral energy 
distribution decreases. The 13 and 33 /im fluxes are synthetic photometry points derived 
from the Spitzer low-resolution spectra, the 70 /im fluxes are photometric data from the 
MIPS instrument on-board the Spitzer Space Telescope. Note that the error bars in the 
upper panel are the same size as the symbols and reflect the small internall uncertainties in 
our IRS spectra. The larger error bars in the lower panel mainly reflect the uncertainties in 
the absolute flux cahbration between the MIPS and IRS instrument. 
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Fig. 7. — The emission bands of crystalline silicates in the 10 /im spectral window. Shown 
are the Spitzer low-resolution spectra normalized to the fitted amorphous silicate, PAH, and 
continuum model (see also Section [2^31) . The normalized spectra all have the same vertical 
scale between 0.97 and 1.15 but are offset for clarity. The spectra are ordered from top to 
bottom by decreasing 10 /im silicate band strength. The ID numbers correspond to those 
listed in Table [H Also indicated in this figure are the positions of the main emission bands 
of silica, forsterite, enstatite dust grains. 
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Fig. 8. — Continuum normalized spectra of the FEPS TTS sample. Shown are the Spitzer 
low-resolution spectra between 17 and 36 /im normalized to the fitted continuum using a low 
order polynomial (see also Section [273]) . The normalized spectra all have the same vertical 
scale between 0.97 and 1.15 but are offset for clarity. The spectra are ordered from top to 
bottom by decreasing 10 /im silicate band strength. The ID numbers correspond to those 
listed in Table [H Also indicated in this figure with the dashed and dotted lines are the 
positions of the main spectral features of forsterite and silica. 
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Fig. 9. — Comparison between the forsterite band position at 27.8 /im and that of the 
33.5 fim feature (top panel) and the 23.4 fim feature (lower panel). Also indicated are the 
calculated nominal band positions for 0.1 fim forsterite grains (filled triangles). 
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Fig. 10. — Comparison between the measured strength of the forsterite bands at 11.3 ^m 
(plotted on x-axis) and 10.1 /im (filled circles), between the 33.5 /xm (plotted on x-axis) and 
the 27.8 /im spectral features (open circles) and the 23.4 //m (plotted on x-axis) and the 
27.8 iim spectral features (filled triangles). 
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Fig. 11. — Correlations between the peak over continuum ratio of the 10 //m sihcate band 
and the band strengths (filled symbols; left axis) and positions (open symbols; right axis) 
of three spectral emission bands of crystalline silicates. The top, middle and bottom panels 
show the correlations for the 23.4 fim, 27.8 fim, and the 33.5 band, respectively. As one 
can see, both position and normalized band strength of these crystalline silicate bands show 
almost no variations and are not correlated to the silicate emission at 10 /im. 
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Fig. 12. — Correlations between the measured band strengths of the crystaUine silicate 
features in the 20 to 35 //m wavelength region and those observed in the 10 jim. spectral 
region. The top figure shows the comparison between the band strength of the 11.3 /im and 
that of the 33.5 /im spectral feature (open circles), the 27.8 /im band (filled circles) and 
23.4 fim band (open squares), respectively. The lower figure shows a similar comparison but 
with the band strength of the 10.1 /im spectral feature. 
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Fig. 13. — Correlation between the enstatite mass fraction of the crystaUine sihcates as 
measured in the 10 /im spectral window and the enstatite mass fraction of the crystalline 
silicates derived from analyzing the longer wavelength features between 20 to 35 /im. 
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Fig. 14. — Observed PAH emission bands in the spectra of the FEPS TTS sample. Shown are 
the Spitzer low-resolution spectra normahzed to the fitted silicate model (see also Section 2.2 
and Table [3]). The first 5 spectra clearly show an emission feature at 8.2 /im we assign to 
emission from PAH molecules. Shown at the bottom of the figure is the average spectrum 
of source # 0,1,3 and 4, over-plotted with the PAH template spectrum used in our fitting 
procedure. The dashed lines mark the positions of the main PAH bands. 



